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ABSTRACT 


We  have  studied  the  SRS  in  liquids.  The  production  of  high  bright¬ 
ness  Raman  devices  is  possible  in  liquids;  however,  the  brightness 
gains  over  existing  laser  devices  does  not  appear  to  be  very  large  We 
must  exercise  care  in  selection  of  the  material  to  achieve  any  semblance 
of  gain. 


The  case  for  gases  is  much  better.  Many  unexplained  processes 
occur  in  H?  and  possibly  other  gases  as  well.  The  brightness  gains 
that  could  Be  produced  have  not  yet  been  determined,  and  the  main 
effort  is  now  concentrated  on  this  problem.  We  are  also  investigating 
the  strange  behavior  of  in  the  usual  SRS  experiments. 
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I. 


INTRODUCTION 


In  principle,  stimulated  Raman  scattering  should  produce  higher 
brightness  sources  than  can  be  obtained  from  conventional  giant  pulfe 
lasers.  For  such  high  brightness  performance,  the  Raman  medium 
must  rernam  essentiaiiv  diffraction  limited  during  the  pumping  time 
in  addition,  the  stimulated  Raman  process  must  proceed  in  a  uniform 
manner  without  transverse  amplitude  or  phase  variations  which  would 
cause  a  reduction  in  brightness.  Initial  experiments  with  liquids  were 
unsuccessful  because  of  the  essential  role  played  by  light  trapping,  and 
as  a  result,  alternative  approaches  were  proposed. 

At  the  beginning  of  the  report  period  we  started  on  a  two-phasi 
program.  The  first  phase  was  to  determine  whether  any  liquids  could 
be  successful  in  the  search  for  a  brightness  gain.  The  second  phase 
involved  a  study  of  the  properties  of  gaseous  systems,  such  as  H-,  and 
CH4,  to  determine  their  characteristics.  The  first  phase  is  now  com¬ 
plete  and  it  appears  that  at  least  one  material,  SiBr4,  will  show  a  bright - 
LeSSegarnu  althou8h  lts  limitations  in  this  respect  still  are  questioned 
Much  of  the  work  on  the  second  phase  is  still  to  be  carried  out,  although 
we  report  here  some  of  our  preliminary  work. 

In  general,  the  details  of  the  work  reported  in  the  earlier  quar- 
erly  reports  are  not  included  here,  but  the  results  are  employed  in 
this  discussion* 
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II  LIQUIDS 


A .  Nonlinear  Index  Effects  in  Liquids 

During  the  past  six  months  we  have  considered  in  detail  the 
nonlinear  effect.*-,  especially  self-focusing,  which  modify  stimulated 
Raman  scattering  in  symmetric -molecule  liquids.  These  liquids 
offered  the  best  hope  for  highly  coherent  Raman  generation.  We  have 
found  that  these  modifying  effects  will  never  be  small,  so  that  for 
example,  self- focusing  will  occur  in  a  symmetric -molecule  liquid  at 
a  threshold  power  not  much  above  that  for  an  asymmetric -molecule 
liquid  having  the  same  particle  density  and  average  polarizability. 

This  is  a  result  of  the  molecular  redistribution  effect,  discussed  in 
detail  in  Appendix  1  (a  reprint  of  an  article  prepared  in  the  summer  of 
1966). 


We  summarize  below  the  developments  since  last  summer  in 
our  understanding  of  the  nonlinear  effects  in  these  reduced  self  focusing 
(RSF)  liquids. 

It  is  stated  in  Appendix  I  that  the  redistribution  (which  uses  the 
Kirkwood  superposition  approximation  of  unknown  validity}  appears  to 
predict  a  nonlinear  index  that  is  an  order  of  magnitude  higher  than 
one  would  guess  from  the  few  measured  thresholds  for  self  focusing. 

For  example  the  theory,  suggests  that  the  symmetric-molecule 
liquid  SiBm  has  a  threshold  almost  as  low  as  the  very  asymmetric 
molecule  CS?.  Attempts  to  measure  this  nonlinear  index  directly  have 
been  initiated  here  and  elsewhere,  but  no  definitive  results  have  yet 
resulted  from  these  attempts.  However,  a  previous  observation  bjr 
J.  L.  Emmett  that  ~5  p  filaments  form  in  stimulated  beams  at  90  to 
the  incident  beam  and  that  these  filaments  have  anomalously  low 
powers,  has  been  reconfirmed.  Perhaps  the  most  significant  measure 
ments  of  nonlinear  indices  to  come  to  our  attention  are  those  of  J. 
Minard,  who  has  measured  directly  the  dc  nonlinear  indices  of  CCl^ 
CS2J  and  benzene  in  each  instance,  he  found  more  index  than  could  be 
accounted  for  by  the  predicted  electrostriction  plus  two  thirds  of  the 
observed  induced  birefringence  (or  Kerr  index1.  In  ea>.h  case,  the 
excess  index  was  around  one  tenth  of  that  expected  on  the  basis  of  the 
estimates  in  Appendix  I  This  finding  is  in  agreement  with  the  non¬ 
linear  index  estimates  from  self  focusing  thresholds  of  an  incident 
laser  beam.  Therefore,  we  now  tend  even  more  strongly  toward  the 
view  that  the  three-  and  four -particle  correlation  functions  are  not 
described  well  for  liquids  by  the  Kirkwood  superposition  approximation 
used  in  Appendix  I.  Just  how  they  are  to  be  described  is  a  problem 
which  we  are  studying. 
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The  saturation  of  the  nonlinear  index  was  not  discussed  in 
Appendix  I,  but  it  is  important  in  determining  for  example,  the  mini 
mum  radius  and  stability  of  a  filament  of  self- focused  light  From 
the  studies  of  the  laser  beam  profile  after  it.  has  passed  through  a  10  cm 
cell  of  SiBr^  and  produced  SRS,  we  conclude  either  that  Raman  scatter¬ 
ing  is  occurring  without  self-focusing,  or  that  the  index  change  saturates 
at  a  value  that  produces  filaments  no  smaller  than  40  to  50  p  We  are 
currently  setting  up  an  experiment  to  measure  the  Raman  scattering 
cross  sections  for  these  liquids  to  determine  whether  SRS  threshold 
could  be  attained  in  our  experiments  without  self-focusing  This  may 
then  determine  indirectly  the  saturation  region 

In  order  to  see  how  assumptions  about  the  short  .range  order, 
other  than  those  tried  in  Appendix  I,  afrect  nonlinear  indices  we  have 
made  calculations  which  assume  that  locally  the  liquid  is  like  a  solid 
with  the  molecules  executing  thermal  and  quantum  vibrations  about  fixed 
sites.  These  calculations  are  not  complete,  but  indicate  that  a  lower 
nonlinear  index  is  expected  in  such  circumstances.  The  nonlinear  index 
is  so  sensitive  to  the  short  range  order  that  it  appears  to  be  a  very  good 
experimental  parameter  for  checking  hypotheses  about  liquid  structure. 


Nonlinear  Experiments  with  Symmetric  Molecule  Liquids 


We  have  conducted  several  experiments  to  study  the  SRS  effect 
in  liquids  composed  of  spherically  symmetric  molecules  Careful 
measurements  detected  no  evidence  of  self  focusing  in  SiBm.  However, 
with  a  multimode  pump,  we  have  observed  evidence  of  a  nonlinear  index 
in  the  form  of  frequency  spreading  in  the  Stokes  radiation.  Threshold 
measurements  for  symmetric  and  nonsymmetric  liquids,  combined 
with  cross  section  measurements,  have  given  a  basis  for  comparison 
of  theoretically  predicted  and  observed  gains. 


1.  Search  for  Self  -  Focusing  in  S^Br^ 

Our  initial  observations  of  the  laser  and  Raman  radiation 
exiting  from  a  10  cm  path  length  of  SiBr^  indicated  some  filamentary 
structure.  However  careful  control  experiments  showed  that  this  was 
probably  a  result  of  a  combination  of  weak  higher  order  transverse  modes 
in  the  GPL,  and  experimental  difficulties  encountered  in  imaging  through 
the  2-m  spectrograph.  In  order  to  obtain  more  conclusive  results  v/e 
set  up  the  apparatus  to  analyze  the  polarization  of  laser  and  Raman  radia  ¬ 
tion,  using  circularly  and  linearly  polarized  pumps  J  This  technique  has 
proved  to  be  a  very  sensitive  test  for  self  •  focusing  in  liquids  like  nitro¬ 
benzene.^  The  radiation  exiting  from  the  cell  is  analyzed  for  polariza 
tion  and  wavelength  and  a  pronounced  depolarization  of  the  radiation  is 
observedto  accompany  the  onset  of  self  -  focusing  (SRS  threshold)  with  a 
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circularly  polarized  pump.  We  detected  no  depolarization  in  SiBr^ 
with  this  apparatus,  although  the  theory  of  self  •  focusing  regardless 
of  the  underlying  mechanism,  indicates  that  there  should  be  strong 
depolarization  if  self-focusing  occurs.  This  negative  result  held  even 
far  enough  above  threshold  to  achieve  substantial  conversion  to  Stokes 
power  and  to  excite  at  least  five  orders  of  Stokes  radiation  It  should 
be  noted  that  by  using  a  1/Z  mm  diameter  pump  by  beam  the  Stokes 
radiation  often  appeared  in  a  small  filament  about  40  to  50  p  in  diam¬ 
eter.  Presumably  this  is  due  to  the  effect  of  spatial  narrowing  oi  the 
Stokes  caused  by  the  variation  of  intensity  across  the  pump  beam  and 
it  attests  to  the  large  gain",  achieved  in  this  liquid. 

2 .  Spectral  Observations 


Using  the  2-m  spectrograph  we  measuied  a  Stokes 
shift  of  225  cm  *  in  SiBr^,  giving  a  first  Stokes  wavelength  of  705  3  mp . 
This  would  be  a  symmetric,  or  "breathing"  mode  of  vibration  of  the 
SiBr^  molecule.  With  tetrabutyl  tin  (TBT),  which  should  also  show 
spherical  symmetry,  we  observed  a  shift  of  2920  cm"*  corresponding 
to  a  Stokes  wavelength  of  871.1  mp .  This  indicates  that  a  C  H  vibration 
has  the  lowest  threshold  in  TBT.  Because  the  higher  order  Stokes 
shifts  in  TBT  fall  outside  the  range  of  the  instrument  only  the  first 


order  Stokes  radiation  was  observed  in  this  material.  With  SiBr^,  up 
to  six  orders  of  Stokes  radiation  were  observed 


Our  giant  pulse  laser  typically  operates  in  several  longitudinal 
modes.  A  Fabry-Perot  interferogram  typically  shows  two  to  four  modes, 
with  spacings  of  about  0.6  cm'*  .  Other  observations  with  the  high 
dispersion  grating  of  the  2-m  spectrograph  indicate,  typically  two 
modes  spaced  by  0.6  to  0  7  cm"  *  .  There  is  usua  ,ly  a  time  modulation 
of  the  pulse,  indicating  other,  more  closely  spaced  modes  probably 
adjacent  longitudinal  modes  of  the  laser  cavity. 

It  has  been  well  established  that  frequency  broadening  is  observed 
in  the  Raman  emission  of  liquids  which  self  focus  when  a  multimode 
pump  is  usedA'  ^  This  is  due  to  a  modulation  of  the  nonlinear  index  at 
harmonics  of  the  difference  frequencies  of  the  pump  ^  We  have  observed 
such  broadening  effects  in  SiBr^  and  TBT,  with  the  amount  of  broadening 
increasing  with  increasing  pump  power  above  threshold,  and  also  increas¬ 
ing  with  the  Stokes  order  for  SiBr^.  We  have  observed  up  to  8.5  mp 
broadening  in  the  fifth  order  Stokes  radiation  from  SiBr^  and  up  to 
7.4  mp  broadening  in  the  first  Stokes  radiation  from  TBT  This 
spreading  of  about  130  cm'*  is  many  times  the  laser  mode  spacing  of 
about  0.7  cm'*,  and  indicates  that  there  is  modulation  of  a  nonlinear 
index  in  these  materials.  However,  at  least  in  Sibr4  this  nonlinear 
index  does  not  lead  to  depolarization  of  a  circularly  polarized  pump,  our 
best  test  for  self-focusing. 
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3. 


Threshold  Measurement s 

Another  good  indioatior  of  whether  self  focusing  is  present 
is  the  theoretical  comparison  of  predicted  gains  using  the  incident  laser 
intensity  with  the  minimum  gain  needed  experimentally  to  observe  Stokes 
radiation.  For  this  purpose,  we  have  made  threshold  measurements  on 
several  liquids.  In  Section  II-C,  below  these  will  be  compar  ed  with 
measured  spontaneous  Raman  cross- sections .  For  the  threshold  mea¬ 
surements,  the  GPL  beam  was  incident  on  a  1.5  mm  aperture.  Follow¬ 
ing  this,  a  40  cm  focal  length  lens  weakly  focused  the  pump  beam  in  the 
10  cm  cell.  A  magnesium  oxide  screen  behind  the  cell  scattered  the 
output  radiation  to  calibrated  Stokes  and  laser  detectors.  The  Stokes 
detection  system  was  capable  of  detecting  about  1  W  of  Stokes  power 
from  the  cell,  except  for  SiBr^  where  the  minimum  detectable  power 
was  somewhat  higher  due  to  the  filters  used.  With  this  arrangement 
the  entire  cell  length  was  within  the  depth  of  focus  of  the  40  cm  lens  so 
that  the  resulting  beam  (about  0.5  mm  ir  diameter)  was  approximately  a 
plane  wave  throughout  the  cell.  The  thresholds  measured  with  this  appara¬ 
tus  are  given  below.  The  error  limits  indicate  the  scatter  in  the  data,  and 
do  not  include  calibration  and  other  systematic  errors.  However  these 
errors,  which  should  be  less  than  20%  should  not  affect  the  relative 
thresholds.  A  linearly  polarized  pump  beam  was  used  for  all  of  these 
measurements . 


Materi, 

Silicon  tetrabromide 
Tetrabutyl  tin 
Nitrobenzene 
Benzene 

Carbon  disulfide 
Carbon  tetrachloride 


Threshold  Power  MW 
0.29  ±  0.02 
0.29  ±  0.03 
0.050  ±  0.005 
0. 14  ±  0.01 
0.024  ±  0.002 
0.65  ±  0.03 


These  powers  should  be  multiplied  by  500  to  give  the  threshol'  intensity 
in  MW/cm^.  The  above  results  for  nitrobenzene  and  benzene  can  be 
compared  to  Wang  s  results-r  by  taking  into  account  the  different  beam 
sizes  and  extrapolating  his  results  to  a  10  cm  cell  length.  We  must 
also  take  into  account  the  fact  that  our  beam  is  nongaussian  in  profile. 

We  have  measured  a  spot  size  in  the  focal  plane  of  the  40  cm  lens 
which  is  four  times  the  diffraction  limit  for  the  1  5  mm  aperture  used. 
Therefore,  if  we  take  a  "characteristic  transverse  radius  of  curvature"^ 
of  the  laser  intensity  of  one  quarter  of  the  beam  radius,  we  find  a  pre¬ 
dicted  threshold  of  0.045  MW  for  nitrobenzene  and  0.16  MW  for  benzene. 
Thus,  our  results  for  these  two  materials  are  in  good  agreement  with 
Wang's  results  and  the  theory  of  self-focusing. 
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C  Cross  Section  Measurements 

In  this  section  we  describe  the  peak  normal  Raman  cross  sec¬ 
tion  measurements  and  compare  the  insulting  theoretical  gains  with 
the  above  threshold  measurements. 

1  .  Discussion  of  Experiments  and  Data 

The  normal  Raman  scattering  data  were  taken  with  a  cw 
argon  laser  operating  at  514.5  mp  and  a  double  monochromator  arrange¬ 
ment  with  photodetection .  The  relative  cross  sections  oi  several  liquids 
were  measured  and  the  absolute  peak  cross  sections  calculated  using 
our  previous  result  of  1  .  3  cm'^  cross  section  per  unit  volume  per  unit 
wavelength  for  nitrobenzene  .6  The  results  are  tabulated  below. 


Liquid 

Raman  l^ine , 
cm- 

Peak  Cross  Section 
? 

cm  “ 

Gain 

cm '  *  /MW/cm1 

Nitrobenzene 

1345 

i  .  3 

1  4  x  10' 3 

Benzene 

992 

2.4 

2.52 

Carbon  Tetrachloride 

440 

0.82 

0.82 

Tin  Tetrachloride 

370 

2.9 

2.54 

Silicon  Tetrabromide 

225 

5.8 

4.23 

Unfortunately,  no  spontaneous  Raman  lines  were  observed  in  TBT  with 
this  apparatus . 

The  accuracy  of  these  measurements  is  not  well  known,  but  ic 
judged  to  be  about  30%.  The  major  limiting  factor  is  that  the  double 
monochromator  used  had  an  instrurr  sntal  width  of  about  0.07  mp  and 
made  it  difficult  to  determine  linewidths  accurately.  A  better  mono¬ 
chromator  will  soon  be  available  which  will  allow  more  reliable 
measurements . 
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Comparison  of  Theoretical  Gain  with  Threshold 

Measurements  i 


a 


>5 


1 

Jm 


The  largest  initial  intensity  available  for  amplify  ation 
in  the  Raman  cell  is  due  to  the  quantum  mechanical  zero-point  vibra¬ 
tions  in  the  electromagnetic  field  and  is  of  order  10‘  ^  W/cm^.  For  our 
beam,  this  represents  2  x  10”®  W  power.  Since  our  detection  system 
was  capable  of  detecting  about  1  W  of  output  power  the  gain  required  is 
about  5x10^.  For  a  double  pass  through  a  10  cm  cell,  this  implies  a 
gain  of  0.88  cm  *  .  Below  we  tabulate  the  theoretical  gains  obtained  by 
multiplying  the  above  gains  (cm'VMW/cm^)  by  the  threshold  laser 
intensity  observed  experimentally. 


! 

i 


r 


Liquid 

Nit  -obenzene 
Benzene 

Carbon  tetrachloride 
Silicon  tetrabromide 


Theoretical  Gair 
0.035  cm  * 

0. 176 
0.266 
0.615 


F 

25 

5 

3.3 

1.4 


F  is  the  factor  required  to  increase  the  laser  power  to  give  a  theoreti¬ 
cal  gain  equal  to  the  gain  experimentally  required  to  observed  Stokes 

emission. 


With  the  many  cumulative  errors  involved  this  comparison 
the  results  appear  adequate  to  explain  the  ODserved  SRS  from  SiBr4  o 
he  tali,  of  purely  exponential  gain.  The  results  for  nitrobenzene  are 
consistent  with  o/r  earlier  measurement..*  and  reaffirm  the  nec...tty 
of  invoking  the  self-focusing  process  to  explain  the  observeo  SRS.  The 
pr«.«  ro'zoK.  for  benzeno'and  carbon  tetrachloride  are  not  conclude 
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III. 


HYDROGEN  GAS  STUDIES 


Preliminary  experiments  were  started  on  H2  gas  to  evaluate  it 
as  an  active  medium  for  a  Raman  laser.  A  50  cm  cell  was  operated 
at  pressures  of  200  to  1000  psig  with  no  detectable  Stokes  output  when 
an  unfocused  GP  laser  was  used  as  a  pump.  Focusing  in  the  cell  with 
a  40  cm  focal  length  lens  produced  an  increase  in  the  power  density 
inside  the  cell  (reported  in  the  section  on  liquids)  and  strong  SRS  out¬ 
put  was  observed.  The  unfocused  output  power  density  of  the  laser  was 
expected  to  be  in  excess  of  40  MW/cm^;  on  the  basis  of  published  data, 
we  felt  that  the  laser  should  have  been  able  to  reach  SRS  threshold  in 
H2  with  no  difficulty.  Tests  with  a  90  cm  cell  did  give  SRS  as  expected 


A .  Near  Field  Power  Distribution  of  the  Laser  Pump 

The  actual  power  density  of  the  laser  beam  at  the  location  of  the 
Raman  cell  was  checked  first.  A  0.  17  mm  diameter  aperture  was 
carefully  scanned  across  the  laser  beam  in  steps  of  0.20  mm  and  the 
laser  output  power  was  compared  with  the  power  measured  on  the 
other  side  of  the  aperture.  The  transverse  mode  structure  of  the  laser 
was  monitored  for  each  data  point.  The  scan  was  made  twice,  once 
with  moderately  successful  efforts  at  transverse  mode  control  and  once 
with  no  effort  to  control  the  transverse  mode.  The  scatter  in  the  data 
was  quite  severe  for  the  non-r.  ode-selected  case  with  noticeable  modu¬ 
lation  on  the  transmitted  beam.  The  mode- selected  data  showed  some 
scatter  caused  by  variations  in  the  effectiveness  of  the  transverse  mode 
control  and  a  much  smaller  residual  scatter.  Systematic  variations  and 
variations  which  are  not  understood  appeared  in  both  cases,  and  again, 
were  much  less  apparent  in  the  mode-selected  case.  The  full  width  at 
half  maximum  for  the  two  cases  (from  the  smoothed  curves)  were  1 .3  ± 
0.2  mm  and  2.9  ±  0.5  mm.  Both  curves  looked  gaussian  in  shape, 
although  some  asymmetries  were  evident.  The  relevant  power  density 
values  were  12  ±  1  MW/crn^*  per  MW  of  power  output  and  3  ±  1  MW/’cm^ 
per  MW  of  power  output.  Typical  power  output  values  were  5  MW 
(mode-selected)  and  10  MW  (non-mode-selected).  These  then  cor¬ 
respond  to  a  power  density  of  60  MW/cm^  (mode-selected  and  30  MW/ 
enm  (non-mode- selected) . 

cm^  (non-mode- selected) . 


B.  SRS  Measurements  in  H2’ 


The  Hi  cell  was  extended  to  90  cm  length  and  the  experiments 
repeated.  SRS  was  observed  at  about  40  MW/cm^  and,  as  had  been 
reported  by  many  other  observers,  it  appears  to  have  a  sharp  threshold. 
In  comparison  with  the  liquids  the  gain  should  have  been  1.5  x  10"^  PQf 
1 . 5  x  40  x  100  x  10  3  6  corresponding  to  the  rather  low  gain  o!  cxp(6). 

The  single  pass  exponential  amplification  from  noise  to  the  1  W  level 
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oi  detectability  of  the  experimental  arrangement  requires  a  gain  of 
exp[17.7]  .  This  is  about  a  factor  of  three  Vo  low  a  laser  power  to 
see  this  signal  and  the  threshold  behavior  suggests  that  there  is  feed¬ 
back  as  has  also  been  reported  by  many  others.  We  have  been  able 
to  observe  some  effects  which  we  think  could  be  feedback  from  the 
end  windows  of  the  ce'1.  Preliminary  data  indicate  that  we  cannot 
consistently  reach  threshold  at  60  MW/cm2  if  the  windows  of  the  cell 
are  tilted  by  as  little  as  0.005  rad  away  from  perpendicular  to  the 
laser  beam.  This  strongly  suggests  that  the  windows  are  contribut¬ 
ing  feedback.  To  prevent  any  feedback  from  the  laser^  and  therefore 
isolate  to  the  single  pass  conditions,  we  installed  a  60  flint  glass 
prism  between  the  laser  and  the  H2  cell.  Other  optics  folded  the 
beam  around  to  again  have  the  laser  aligned  into  the  H2  cell  but  now 
no  feedback  to  the  laser  could  occur  at  the  Stokes  wavelength  due  to 
the  dispersion  in  the  prism.  To  our  surprise  the  threshold  behavicr 
still  occurred  -  ai  about  40  MW/cm2.  This  result  has  not  ruled  cut 
feedback  from  the  windows  of  the  Raman  cell  but  the  prism  has  pre¬ 
vented  feedback  from  the  laser  cavity.  Also  it  was  known  that  the 
laser  beam  was  not  perpendicular  to  either  window  of  the  H2  cell  when 
the  dispersing  prism  was  installed.  Further  tests  of  the  above  resu.us 
are  contemplated,  although  all  apparent  sources  of  feedback  have  been 
removed. 
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IV. 


CONCLUSIONS  AND  PLANS  FOR  THE  FUTURE 


The  most  promising  approach  to  achieving  high  brightness 
appears  to  be  the  use  of  a  Raman  material  which  will  retain  good 
optical  properties  during  the  stimulated  scattering  process.  Although 
satisfactory  liquids  may  be  found  eventually,  the  gases  hold  the 
greatest  promise  of  brightness  gains.  We  are  currently  installing  a 
large  ruby  amplifier  in  order  to  pump  a  Raman  oscillator -amplifier 
combination.  The  Raman  amplifier  cell  is  under  construction  now 
and  will  have  a  clear  aperture  of  l-l/2  in.  diameter  with  demount¬ 
able  sections  for  a  5,  10,  15,  or  20  ft  length.  We  are  currently 
working  on  the  lab  arrangement  to  accommodate  the  oscillator- 
amplifier  cells.  We  plan  to  make  measurements  of  the  output  beam 
characteristics  from  the  amplifier  under  various  operating  conditions. 
Until  this  apparatus  is  installed  we  will  continue  with  the  measure  • 
ments  on  the  H2  oscillator  in  an  effort  to  achieve  a  better  understand¬ 
ing  of  the  H2  gas  behavior,  which  still  presents  disturbing  anomalies. 

With  the  possible  exception  of  SiBr^,  SRS  in  the  liquids  we  have 
studied  seems  to  offer  little  hope  of  achieving  brightness  gains  over  a 
conventional  GP  laser.  Even  in  SiBr^,  we  expect  that  self -focusing  will 
occur  at  power  levels  only  "omewhat  larger  than  those  reouired  for 
other  liquids,  such  as  benzene.  Moderate  increases  in  tightness  of 
up  to  about  a  factor  of  10  could  be  obtained  in  short  pulses  by  two 
techniques.  First,  if  a  GP  laser  is  mode  locked  the  power  output 
increases  significantly.  If  there  is  no  degradation  of  the  transverse 
mode  structure,  there  will  be  a  corresponding  brightness  increase. 
Second,  one  could  use  the  large  amplitude,  short  duration  backward 
traveling  pulse  in  a  Raman  oscillatoi  7  to  achieve  the  same  effect. 

These  short  pulses  would  be  useful,  for  example,  in  high  power  propa¬ 
gation  studies. 
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Effect  of  Molecular  Redistribution  on  the  Nonlinear 
Refractive  Index  of  Liquids* 

R.  W.  Hkllwabth 

Hughes  Research  Laboratories,  Malibu,  California 
(Received  6  July  1966) 

An  expression  (or  the  static  nonlinear  dielectric  constant  of  r  nonpolar  liquid  is  derived  which  takes  into 
account  the  reorientation  and  local  spatial  redistribution  of  molecules  in  the  presence  of  a  strong  electric 
field.  This  result  is  used  to  calculate  the  nonlinear  index  of  refraction  of  various  lossless  liquids  at  "optical” 
frequencies  which  are  much  higher  than  molecular  reorientation  rates.  From  this  nonlinear  index,  we  obtain 
estimates  of  the  optical  power  required  to  initiate  self-focusing  in  a  liquid  over  periods  so  short  that  macro¬ 
scopic  density  changes  do  not  have  time  to  take  place.  This  situation  is  commonly  approached  in  experi¬ 
ment.  The  theory  predicts  that  the  effects  of  molecular  redistribution,  hitherto  not  considered,  will  be 
generally  as  important  as  the  widely  considered  effects  of  molecular  reorientation,  and  will  in  fact  dominate 
the  nonlinear  index  of  many  symmetric  and  nearly  symmetric  molecules.  Furthermore,  the  results  suggest 
that  some  purely  symmetric  molecules,  such  as  SiBr»,  may  exhibit  self-focusing  in  liquid  as  readily  as  do 
some  commonly  studied  asymmetric  molecules,  such  as  nitrobenzene.  The  calculation  proceeeds  trom  classi¬ 
cal  statistical  mechanics  with  the  aid  of  a  variational  principle  that  is  valid  for  arbitrary  density.  In  the 
low-density  limit  our  results  reflect  only  molecular  reorientation  and  reduce  to  those  of  Debye  and  others. 
The  accuracy  of  the  results  depends  mainly  on  the  accuracy  of  Kirkwood’s  “superposition  approximation” 
in  representing  three-  and  four-particle  correlation  functions  in  liquids.  Since  the  accuracy  of  this  approxi¬ 
mation  is  at  present  unknown,  nonlinear  index  data  may  prove  useful  in  checking  it.  As  a  by-product  of  our 
investigation,  we  Lave  proven  that  the  Clausius-Mossotti  expression  gives  a  lower  bound  for  the  dielectric 
constant  of  a  fluid  in  which  the  two-particle  correlation  function  is  a  function  only  of  the  interparticle  spacing 
and  approaches  a  constant  at  large  spacing.  Molecular  redistribution  must  also  play  a  role  in  induced  bire¬ 
fringence  (ac  and  dc  Kerr  effects),  especially  for  symmetric  molecules.  However,  the  present  treatment  is 
limited  to  waves  of  a  single  linear  polarization,  and  does  not  cover  the  Kerr  effect. 


I.  INTRODUCTION 

THE  self-focusing  of  optical  beams  and  the  effects 
of  this  self-focusing  on  stimulated  Raman  scat¬ 
tering  and  other  nonlinear  effects  have  been  studied 
in  a  variety  of  liquids.1-8  The  nonlinear  index  of  re- 

•  This  work  has  been  supported  in  part  as  a  part  of  Project 
defender  under  the  joint  sponsorship  of  the  Advanced  Research 
Projects  Agency,  the  Office  of  Naval  Research,  and  the  Depart¬ 
ment  of  Defense. 

1 N.  F.  Pilipetskii  end  A.  R.  Rustamov,  JETP  Pis’ma 
Redaktsiyu  2,  88  (1965)  [English  transl.:  JETP  Letters  2,  55 
(1965)]. 


fraction  which  causes  self-focusing  has  been  attributed 
to  electrostrictive  effects  and  to  the  reorientation  of 
asymmetric  molecules  (ac  Kerr  effect)  in  the  strong 
optical  fields;  the  electronic  nonlinearities  are  presumed 


*  G.  Hauchecome  and  G.  Mayer,  Compt.  Rend.  261,  4014 
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tu  eu.itnbuto  iinuli  less.-  hi  this  |ki|k t  w.-  point 
,uit  th.it  .mother  meihankm.  the  spatial  redistribution 
of  the  moles  iik"'.  i an  male  i  contribution  to  the  non¬ 
linear  index  of  liquids  whit  ti  is  generally  comparable 
U)j  |  nla\  in  some  eases  <lominale,  eke t rostra  tive 

ami  reoiieniatioil  e'aeets. 

\\e  will  consider  only  the  index  change  experienced 
l,v  a  strong  I i hoi rly  polari/ed  optical  neld  in  a  liquid; 
the  index  changes  for  beams  of  other  polarizations  (and 
related  phenomena  stall  as  induced  birefringence)  wilt 
involve  added  cakulalional  lomplesi'ics.  The  strong 
impressed  held  induces  a  dipole  moment  in  each  mole- 
u,le  so  that  extra  anisotropic  dipole  dipole  forces,  which 
are  evidentk  proportional  to  the  tune  average  of  t  ie 
square  of  the  field  strength,  oeetir  between  molecules 
in  addition  to  the  usual  intermolcetdar  forces  already 
present,  Therefore,  after  the  molecules  come  into 
thermal  equilibrium  in  the  impressed  held,  the  n- 
particle  correlation  functions  for  the  liquid  («>  *)  are 
altered  in  an  anisotropic  way  from  their  zero  held  forms. 
That  is,  the  field  causes  the  molecules  to  become 
“redistributed”  in  space,  and  this  redistribution  con¬ 
tributes  to  a  change  in  index.  For  example,  the  induced 
dipole  forces  cause  two  nearby  molecules  of  a  fixed 
separation  to  be  more  probably  found  in  line  with  the 
strong  applied  field  (where  their  combined  polarization 
is  maximum)  than  on  a  line  perpendicular  to  it  (where 
their  combined  jrolarization  is  minimum).  Of  course, 
asvmmctrie  molecules  also  tend  to  reorient  themselves 
in  the  field,  and  thus  increase  their  average  polartza 
bilitv,  and  the  macroscopic  mass  motion  caused  by 
electrostrictive  forces  and  thermal  expansion  also  alters 
dielectric  properties  after  macroscopic  density  vari¬ 
ations  have  had  time  to  form.  It  is  the  mam  object  of 
this  paper  to  assess  the  change  in  index  of  refraction 
brought  about  bv  the  local  redistribution  and  reorien¬ 
tation  of  molecules  in  a  liquid,  ignoring  the  effects  of 
macroscopic  density  changes  which  take  much  more 
time  to  form. 

Specificallv,  we  calculate  here  the  nonlinear  index 
(to  second  order  in  the  fields)  which  would  exist  in  a 
fluid  irradiated  by  an  optical  pulse  of  such  short  dura¬ 
tion  that  macroscopic  changes  in_ density  (such  as 
might  arise  from  electrostrictive  effects  and  heating) 
do  not  have  time  to  occur.  We  also  assume  that  a 
molecule  will  undergo  so  many  col'isions  during  mis. 
optical  pulse  that  a  statistical  mechanical  equilibrium 
will  exist  which  is  appropriate  to  the  modified  inter- 
molecular  dipole dipole  forces  and  the  torques  resulting 
from  the  application  of  the  strong  optical  fields.  The 
refractive  index  for  this  altered  equilibrium  state  of 
the  fluid  is  then  calculated  in  a  standard  way  by 
differentiating  the  free  energy  a  of  a  uniform  sample 
of  fixed  macroscopic  density  with  respect  to  the  applied 
electric  field  strength  to  obtain  the  induced  (hpole 
moment.  From  this  one  obtains,  in  the  usual  way,  --.n 
expression  for  the  total  (linear  plus  nonlinear)  index. 

~Tjp s'-.  Chino,  1-  Garmke.  and  C.  H  Townes,  Phys.  Rev. 
tellers  13,  479  (1961). 


This  expression  for  the  index  contains  only  the  effects 
o.'  molecular  redistribution  and  reorientation. 

The  densities  of  the  liquids  of  interest  are  too  large 
for  a  low  densitv  expansion  of  d  to  be  accurate.  There 
fore  we  employ'  a  variational  principle  which  gives  an 
upper  bound  for  5  regardless  of  density.  A  local  field 
approximation”  is  then  made  in  which  certain  free 
parameiers  may  be  varied  to  obtain  a  lowest  upper 
bound  for  5  within  the  local  field  approximation.  From 
this  a  lower  bound  results  for  the  linear  part  of  the 
index  of  refraction,  which  bound  is  the  “Clausius- 
Mossotti”  or  “Loren tz  Lorenz”  expression.  This  is  a 
strong  lower  bound,  true  whenever  two  particle  corre¬ 
lations  depend  only  on  the  interparticle  spacing.  Ln- 
fortunately,  a  similar  bound  does  not  result  for  the 
nonlinear  index  which  we  calculate  to  second  order  in 

the  field.  ,  ,  , 

Because  the  Clausius-Mossotti  expression  for  die 
linear  dielectric  constant  (or  the  square  of  the  index) 
is  known  to  be  accurate  to  within  a  few  percent  for 
liquids,  one  might  expect  our  corresponding  expression 
for  the  second  order  part  of  the  index  to  be  accurate  to 
within  10  to  50%.  Unfortunately,  many  additional 
errors  are  introduced  into  the  evaluation  of  the  non¬ 
linear  index  expressions  by  our  having  to  use  tie 
Kirkwood  “superposition  approximation  for_,U'° 
three  and  four-particlc  correlation  functions.  Ihese 
errors  are  of  uncertain  magnitude  within  our  present 
knowledge  of  statistical  functions  for  liquids.  Perhaps 
the  nonlinear  index  experiments  can  be  used  to  obtain 
an  dea  of  the  accuracy  of  the  Kirkwood  approximation, 
sir  ic  the  nonlinear  index  is  more  sensitive  to  the  char¬ 
acter  of  third-  and  fourth  order  correlation  functions 
than  are  other  quantities  which  are  as  accurately  meas¬ 
urable.  It  is  reassuring  that  the  term  in  our  final  ex¬ 
pression  for  the  nonlinear  index  which  dominates  at 
low  densities  and  which  is  most  sensitive  to  the  asym- 
metry  of  the  molecules  is  exactly  the  same  as  the 
expression  of  Debye9  for  the  low-frequency  nonlinear 
index  of  molecules  having  no  permanent  dipole  moments 
(these  don’t  interact  with  optical  fields).  However, 
optical,  rather  than  lov. -frequency,  polarizabilities 
appear.  Although  this  “Debye”  term  m  the  nonlinear 
index  is  dominant  (for  asymmetric  molecules)  at 
densities  somewhat  lower  than  liquid  densities,  the  other 
terms  we  find  arc  important,  if  not  dominant,  at  liquid 
densities.  These  other  terms  arise  from  the  redistri¬ 
bution  of  molecules  brought  about  by  the  field-induced 
changes  in  intermolecular  forces. 

The  relative  power  thresholds  for  producing  seD 
focusing  or  “self  trapping”  in  those  liquids  for  which 


•  P.  Dcbve,  Marx'S  llandbiuh  der  Radiolasie  17  {Leipzig.  1025;, 
Chan  V  n  76S.  P  Ungcvin,  J.  J  Urnv>r.  M.  Born  and  others 
carder  derived  similar  rcla.i-ns  i„  .hat  of  Debye,  but  none  had 
evaluated*1  correctly  both  ihe  local  fields  P*,r«  non- 
linearities  and  lhc  ciTrd  of  Ihe  change  in  .olari/abil  l j  ot  a 
molecule  on  the  local  fields  ihcmselvcs.  furthermore,  in  OMi * 
by  M  Born  (Springer-Yerlag  Berlin  1%5>1  “  ' oTti  366 
is  inln'Uucct!  in  passinj:  from  I  «|-  W),  P-  ^  2  1  Indices 

which  is  copied  in  some  of  ihe  later  literature  on  no'ihnear  indices 
and  lhc  Kerr  effect.  See  NI.  Bom, » bid.  for  a  general  bibliography. 
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quantitative  lata  1^  avail, till*-  arc  consistent  with  values 
cakulaicd  from  our  nonlinear  indites  and  l lie  formula 
for  tin; shold  [lower  of  Chian  ct  <il.‘  In  paititular,  the 
recent  observation  by  Brel  ct  ill."  that  the  symmetric 
molecule  COi  exhibits  a  threshold  anomaly  character¬ 
istic  of  self-focusing  in  liquid  when  the  asymmetric 
molecules  eholofonn,  acetone,  and  acetic  acid  do  not  is 
consistent  with  our  predictions.  Comparison  of  our 
results  w  ith  the  current  data  indicates  that  our  predicted 
relative  threshold  [lowers  may  lie  accurate  to  around 
±25%.  If  this  is  so,  then  we  predict  that  several  purely 
symmetric-molecule  liquids,  such  as  SiBrr,  will  exhibit 
low  self  trapping  thresholds  of  the  same  order  as,  for 
example,  that  of  the  strongly  self  focusing  liquid  nitro¬ 
benzene.  Molecular  redistribution  must  also  affect  the 
induced  birefringence  (Kerr  effect),  perhaps  accounting 
for  the  general  disagreement  between  the  usual  Debje 
theory*  and  experiment  (especially  pronounced  for 
symmetric  molecules).  However,  we  do  not  consider 
such  effects  here. 

The  absolute  magnitudes  of  our  predicted  thresholds 
tend  to  fall  an  order  of  magnitude  lower  than  those 
estimated  from  experiments.2  6  The  errors  inherent  in 
the  Kirkwood  approximation  might  well  be  responsible 
for  this.  However,  as  we  point  out  later  in  the  discussion, 
there  is  some  evidence  that  thresholds  are  actually 
lower  than  was  at  lirst  supposed.  The  possible  effects 
of  relaxation  processes  which  may  prevent  the  achieve¬ 
ment  of  the  quasi  steady  state  vve  have  assumed  have 
not  yet  been  examined,  nor  has  the  role  of  bubbles, 
impurities,  normal  thermal  turbulence,  etc. 

II.  FORMULATION 

The  induced  orientation  and  intermolecular  potentials 
are  quadratic  in  the  field  strength.  Therefore  molecules, 
which  can  neither  translate  nor  rotate  at  optical  fre¬ 
quencies,  will  assume  the  same  configuration  in  a 
uniform  dc  applied  field  as  in  a  uniform  optical  field  of 
equal  mean  square  amplitude  (if  the  optical  and  static 
polarizabilities  are  equal).  We  will  take  advantage  of 
this  fact  to  simplify  the  mechanics  of  the  calculation 
and  compute  the  nonlinear  dielectric  constant  to  second 
order  in  the  fields  as  if  the  strong  applied  electric  field 
were  static.  Of  course,  only  induced  polarizations  need 
be  considered ;  there  is  essentially  no  interaction  be¬ 
tween  any  permanent  electric  dipole  moments  and  the 
optical  fields  and  polarizations.  The  desired  second 
order  optical  dielectric  constant  is  then  simply  obtained 
from  the  static  constant  by  replacing  the  square  of  the 
static  field  by  the  mean  square  optical  field  and  us.ng 
the  optical  values  for  molecular  polarizabilities  in  all 
expressions. 

To  compute  the  static  dielectric  constant,  we  first 
compute  the  free  energy  J  of  an  ellipsoidal  sample  of 

*  G  Bret  F  Gins,  ami  O.  Mayer,  in  I  I  I  I'  J  Quant.  MUlron. 
QR-2,  No.  4,  18  (I960). 


uniform  fluid  in  a  uniform  x-directcd  electric  field  Fa, 
with  a  principal  axis  of  tile  ellipsoid  parallel  to  the  field. 
Then  the  total  electric  dipole  moment  of  the  sample  M 
is  obtained  from  the  well  known  relation10 

M=-,lJ/JE0.  0) 

Because  the  sample  is  ellipsoidal  and  the  external  field 
uniform,  the  polarization  per  unit  volume  P  inside  the 
sample  will  be  uniform  and  x-directed.10  We  assume 
the-  sample  to  have  a  fixed  volume  V,  whence 

P-M/V.  (2) 

We  desire  the  nonlinear  index  only  to  second  order  in 
the  electric  fields  so  we  compute  to  fourth  order  in 
the  fields  and  obtain  a  result  of  the  form 

-ff=(W<oV2+c/-o4/4)l',  W 

from  which  we  have 

P=bEa+cFa f.  (4) 

In  an  ellipsoidal  sample  uniformly'  polarized  along  a 

principal  axis,  the  relation  of  the  polarization  density 
to  the  macroscopic  electric  field  E  is  generally  written10 

E=E0-LP,  (s) 

where  L  is  called  the  depolarization  factor  for  that  axis 
of  the  ellipsoid.  The  factor  L  can  have  values  between 
0  and  4 it  depending  only  on  the  shape  of  the  ellipsoid.10 
If  we  define  a  linear  susceptibility  X  and  a  nonlinear 
susceptibility  ij  by  the  relation 

P=XE+n&,  (6) 

then  a  comparison  of  (4)  and  (6)  with  the  aid  of  (5) 
shows  that 

X=b/(\-Lb)  (0 

and 

r,=c(\-Lb)-*=c(l+Lxy  (8) 

The  nonlinear  dielectric  constant  t.vt  is 
€,v/,=  l-\-4irP/E 

=  l+4ir(X+»jZ?)  (9) 

and,  in  extrapolating  our  results  to  the  optical  case, 
vve  would  therefore  have  an  index  of  refract. on  t.vc15 
given,  to  second  order  in  E,  by 

t,v  l1,j = «+  (2tr/ tP/yE1 ,  (10) 

where  £  is  the  root  mean  square  electric  field  a  veraged 
over  many  optical  cycles  and  11  is  the  ordinary  linear 
index  of  refraction.  Following  the  standard  practice, 
we  define  a  “nonlinear  index”  ni  by 

2jp}/«  .  (11) 

Then  the  threshold  power  Pc  required  to  begin  trapping 

a  uniform  cylindrical  beam  much  broader  than  a  wave- 

10  W  i'  Brown,  Jr.,  llandbticli  dcr  Plirsik  -VI  / /,  Dielectrics, 
edited  by  S.  I'luggc  (Springcr-Verlag,  Berlin,  1956),  p;>  l-t.x3. 
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length  i*  given  approximately  bv5 

r,-  (].->2v.v„,  (ias//-).  (12) 

where  \0  and  cg  are  llie  v;u  mint  wavelength  and  velocity 
of  light,  respeetix  ely.  Wc  now  have  all  the  relations 
necessary  to  council  the  noulittc.tr  index  effects  with 
the  electric  free  energy  ;T. 

In  computing  the  free  energy  A,  we  will  assume  that 
we  are  dealing  with  a  classical  llttid  and  use  chissical 
statistical  mechanics.  We  must  therefore  formulate  the 
total  potential  energy  ;  of  the  .V  molecules  in  the  sample 
when  they  are  plated  at  positions  r‘,  ■  •  •  rv,  have 
electric  tli|K>lc  moments  m1.  ••  •  m  \  and  arc  oriented 
at  the  Euler  angles  0t,  $,,$■„  •  •  ■ .  8\,  <p.\,  ts,  which  wc 
symbolize  by  lb,  •  •  •,  12.v- 

We  may  anticipate  that  the  electronic  nonlinearities 
will  be  small  compared  to  those  with  which  we  are 
concerned  here  and  assume  that  the  internal  energy  rT 
of  the  7th  molecule  is  quadratic  in  the  components  of 
its  elec  trie  dipole  moment  mT.  If  w.’  arc  the  components 
of  this  moment  along  the  principal  axes  <i=l,  2.  2,  of 
its  polarizability  ellipsoid,  then"1 

n«,  (13) 

a-1 

where  a,,  aj,  aj  are  the  polarizabilities  of  the  molecule 
along  its  principal  axes.  The  components  may  are 
related  to  thecomponents  m,y  measured  in  a  laboratory 
fixed  coordinate  system  by" 

m,y- R,a(Uy)may.  (14) 

The  repeated  space  indices  here  and  elsewhere  are 
assumed  to  be  summed.  The  matrix  R,a(Uy)=R  "tm,(Uy) 
represents  the  rotation  operator  for  the  Euler  angles 
of  the  7th  molecule,  md  will  be  abbreviated  R,ay. 

lire  interaction  energy  vya  between  two  molecules  at 
positions  rT  and  r“  will  be  assumed  to  be10 

r-n.=itv.+m,ym,‘‘D,,y“  (15a) 

where  in-,*  is  a  hant  core  repulsive  potential  and  the 
remaining  term  is  the  dipole  dij>ole  interaction  energy. 


The  may  which  measure  the  displacement  of  electronic 
charge  may  be  taken  to  be  the  independent  harmonic 
“internal”  coordinates  of  the  molecules  which  describe 
the  electronic  polarization.  That  is,  it  is  appropriate  to 
integrate  over  all  values  of  these  coordinates  in  cx’alu- 
ating  the  free  energy  from  the  potential  energy  func 

11  We  will  use  ihe  lower  case  I.atin  subscript  a  ihrough  h  lo 
denote  the  space  components  of  a  vector  or  matrix  referred  jo 
the  principal  axes  (1,2,3)  of  the  polarizability  ellipsoid  fixed  in 
the  molecule  I.atin  subscripts  from »  onward  wifi  be  used  to  denole 
space  com|Hinenls  referred  to  Ihe  orthogonal  principal  axes  (r,v,r) 
of  the  ellipsoidal  sample  fixed  in  the  laboratory. 


lion.10  The  fret-  energy  5  of  the  liquid  in  the  ad  i  retted 
external  field  may  therefore  he  written 

</{r  (16a) 

where 

(16b) 

is  the  total  energy  of  a  configuration  of  the  positions, 
moments,  and  orientations  of  the  molecules,  f  J[i) 
symbolizes  the  integration  fdt1-  •  •fdt*  over  the 
entire  volume  1'  of  the  ellipsoid  of  every  molecular 
position.  Similarly,  </\/{U}  symbolizes  the  integral 
fa’  S\\\dxdd\fa7’  d<t>ifa7’  #|,  •  •  fa'  smdsdd\fa7’ 
Xf>"  #a’(8*5)  ,v  over  all  molecular  orientations  [V.), 
and  ftl[ma)  symbolizes  the  integral  fdmx'dnh'dmj'-  ■  • 
f dmfdmfdmf  over  all  internal  molecular  coordi¬ 
nates  {wa}.  0  is  defined  as  the  inverse  of  (the  tem¬ 
perature  times  Boltzmann’s  constant). 

This  completes  the  formulation  of  the  static  dipolar 
free  energy  suitable  for  classical  liquids.  The  relations 
(1 )  to  (12)  connect  the  dependence  of  5  on  an  externally 
applied  electric  field  E0  with  the  threshold  power  for 
self-focusing  or  self- trapping  which  we  wish  to  predict. 
It  remains  to  carry  through  the  evaluation  of  the  free 
energy’  starting  from  the  basic  definition  (16). 

m.  A  VARIATIONAL  PRINCIPLE  FOR 
CALCULATION 

In  order  to  obtain  from  (16)  an  accurate  estimation 
of  5  for  liquids,  it  does  not  suffice  to  expand  5  in  powers 
of  the  number  density  p= A*/  V.  We  therefore  employ  a 
variational  approach  which  obtains  an  upper  bound 
for  IT.  Parameters  will  exist  in  the  upper  bound  which 
can  be  varied  to  obtain  a  low'est  upper  bound  within 
any  approximation  scheme. 

To  obtain  a  variational  principle,  we  note  that  the 
total  potential  energy  v  is  a  quadratic  function  of  the 
dipole  moment  coordinates  may.  This  means  that  if  we 
define  new  coordinate  variables  of  integration  vay  by 

i'uT=»i«T— ^«T,  0?) 

where  the  n„y  arc  those  values  which  minimize  v  for  a 
given  positional  and  orientational  configuration,  then 
wc  can  rewrite  (16)  as 

r*T=  j  e*i>(-0L».ir*>.“/fc' 

— 0  mill;,:#) ,  (I**) 

where" 

*=*({,.), (r), (!»},/:„) 

~lly(-^iyEn+  \  /Ja15 

n«-fj  ,  (1») 
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is  the  function  whiih  is  minimized  hy  our  choice  of  the 
uay  and 


c\|*[- 


XexpC  3d  °u4" !C»  piy } ' 


(20) 


The  symbol  {y}  represents  the  set  of  all  values  of  the 
PaT  The  intennolecular  potential  t'  defined  by  (20) 
cannot  depend  on  the  electric  field  £<>  (or  the  (y)  would 
not  minimize  <t>)  and  represents,  of  course,  the  long 
range  van  der  Waal’s  potential  due  to  induced  dipole- 
induced  dipole  interactions.  This  potential  has  been 
analyzed  in  the  literature’0  and  we  need  not  discuss  it 
further  here. 

It  is  too  difficult  to  find  an  exact  explicit  expression 
for  [minj^i^],  but  if  we  insert  wrong  values  for  the 
<(•}  into  <p,  we  at  least  obtain  an  estimate  for  the  field 
dependent  part  Sc  of  5  that  is  too  high.  We  therefore 
have  the  useful  inequality  that  for  any  (y)  whatever 


where 

c~BF= j  Jd{T)d{a)uHt),mrB* 

and  „ 

e-\p(-flE„§»0T.-fr') 

fd{r)d{n)  exp(-/3i;T,  Ji\>Y„-/Si') 


(21) 

(22a) 

(22b) 


With  this  inequality  we  may  proceed  with  a  variational 
calculation  of  the  desired  field-dependent  part  SFg  of 
the  free  energy. 


IV.  THE  EFFECTIVE  FIELD  APPROXIMATION 

We  will  use  a  local  field  approximation  in  which  we 
use  the  approximate  values 

(23a) 

or  equivalently" 

(23b) 

in  <t>  as  if  each  molecular  moment  saw  only  an  effective 
x  directed  local  field  of  magnitude  £'.  We  will  then 
adjust  £'  to  minimize  the  resulting  upper  bound  F  for 
the  free  energy  (calc'ated  to  fourth  order  in  £'). 
The  F  thus  obtained  is  a  function  of  the  shape  of  the 
ellipsoidal  sample,  but  the  dielectric  constant  derived 
from  it  via  (9)  is  not.  (Otherwise  one  could  vary  the 
shape  parameters  to  improve  further  the  accuracy'  of 
the  estimate  for  the  dielectric  constant.) 

Using  (23)  in  (19)  gives 


(24) 

iqs^jG'^'E', 

(25a) 

G^Eo— }E', 

(25b) 

1<5=>  £t.,E'  D”  «r  E'. 

(26) 

We  have  introduced  the  usual  shorthand  notation  for 
multiplication  among  the  3X3  space  matrices 
a>s£,aTa./?>.T,  and  the  3-vectors  E0=f£o 

and  E'=i£'.  The  double  average  brack. ts  ((  ))  in 
(24)  symbolize  the  integrations  over  {r}  and  {0} 
weighted  by  ir({r}j  which  is  normalized  so  that 
({!))=  1.  We  will  often  use  only  the  single  average 
bracket  when  cither  the  angle  or  position  average  is 
required  alone;  which  average  is  intended  will  be  clear 
from  the  context.  In  either  case,  <1)—  1. 

We  expect  (as  we  will  find)  that  the  best  local  field 
£'  will  be  of  the  order  of  £0.  Therefore  (24)  implies 
that,  to  fourth  order  in  £0,  F  is  given  by 


F=  Ut—U  i— jj— 2  U  \r\-U  u] 

(27) 

where 

(28a) 

and 

U,i=  ((«,*>)) -  ««.»««/» =  *,  J = 1,2. 

(28b) 

The  integrals  in  t/,  and  U u  are  easily  performed  to 

give 

Vi=NaGE 

(29) 

and 

Uu=  N&E'WA 

(30) 

whtre  the  dimensionless  anisotropy  parameter  A  is 
defined  by 

is 2[(ai—  aj)s+  (aj— aj)’+  (ai—  Oi)*l/ (45a*)  (31) 

and  (a1+aJ-i-ai)/3  is  the  usual  linear  molecular 
polarizability. 

To  perform  the  integrals  in  U >  and  Uj»  we  require 
the  positional  average  of  We  assume  the  weight 

function  w  does  not  depend  on  {0}.  Then,  if  the 
dimensionless  two-particle  correlation  function  p(r) 
defined  by 

p(r»)^V*JdTydrM{r))  (32) 

approaches  a  constant  value  when  rlt=  |n— r,|  is  still 
small  compared  with  the  size  of  the  liquid  sample 
under  consideration,  the  required  positional  average 
gives  (neglecting  terms  of  order  N~l  smaller) 

<Ly..D„T‘)  =  Npl  (33) 

where 

/s£-4t/3,  (34) 

independently  of  the  form  of  p{r).w  p  is  the  number 
density  of  molecules  N/V;  L  is  the  depolarizing  factor 
along  that  principal  axis  of  the  ellipsoidal  sample  which 
is  parallel  to  the  external  field.  Therefore 

l/*=iATo»£V  (35) 

and 

U[i—  .Vo*6'£'VA .  (36) 

To  perform  the  integrals  required  in  l'n,  the  exact 
form  of  the  two-,  three-,  and  four  particle  correlation 
functions  is  required.  For  purposes  of  calculation,  we 
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have  assumed  lh.it  p(r)  is  a  simple  step  hint  tion  whit  h, 
properly  normalized,  is 


where 


p(r)=0,  r<il 

-  (1-r  1')  ',  r>,l 

T  =  4ir(/y.l. 


(’7a) 


(371>) 


We  have  used  the  usual  "Kirkwood  supeqrosition 
approximation”  for  the  three-  and  four  particle  corre¬ 
lation  functions,  expressing  them  in  terms  of  p(r).  The 
details  of  the  calculation  of  l’n  with  these  approxi¬ 
mations  is  given  in  the  Appendix;  the  result  is 

U„/  (,Y£VV)  =  (P+ 2jr*/3)A+  2»»/3 

+  trs(32+136A+167As)/(+V) 

-fO.%lPr(t+AV4),  (38) 


where  A«  is  a  small  contribution  (^lO-1)  from  an 
integral  that  we  have  not  been  able  to  evaluate  exactly 
and  which  we  may  neglect  without  affecting  the 
probable  accuracy  of  the  over-all  results. 


V.  THE  BEST  EFFECTIVE  FIELD;  LINEAR  AND 
NONLINEAR  INDEX  FORMULAS 


Having  explicit  expressions  for  all  the  terms  in  (27) 
for  F,  we  will  adjust  the  magnitude  £'  of  the  effective 
local  field  to  obtain  as  low  a  value  as  possible  for  the 
upper  bound  F  of  the  electric  free  energy.  Since  the 
terms  of  order  £0*  arc  supposed  to  be  orders  of  magni¬ 
tude  smaller  than  those  of  order  £0J,  it  is  easily  seen 
that  the  best  E'  is  that  which  minimizes  the  terms  of 
order  £0J.  Any  small  change  of  E'  front  such  a  value 
raises  these  large  terms  more  than  it  can  depress  the 
smaller  terms  of  order  £04.  That  is,  we  need  only  mini¬ 
mize  (i’t—Ui)  with  respect  to  F'  and  find  that 

£'=£0/(l+po/)  (39) 

is  the  best  choice  for  £'.  Using  (39)  in  (27)  gives  im¬ 
mediately  the  approximate  b  and  c  coefficients  of  (7). 
Since  the  inequality  (21)  holds  in  the  limit  of  low  fields, 
the  estimate  for  b  is  a  lower  bound.  Then,  from  (7)- 
(12),  the  following  relations  result  for  the  linear  di¬ 
electric  constant  e  and  the  nonlinear  index  ns: 


and 


e> 


l+8trpa/3 

1  —  4 rpa/ 3 


(40) 


g-pa’/S  r  4irpa! 

- - A+ - (32+  136A+  167A!) 

n  ( 1  — f  jrpa/3)  t  45  t 

8*Va*(l  +  A)  "I 

+ - +3.843aVr  .  (41) 


We  reiterate  that,  since  the  test  sample  of  fluid  con¬ 
sidered  was  uniformly  irradiated  and  fillet1  a  fixed 
volume  uniformly,  the  estimate  (41)  for  the  nonlinear 


Tahir  i.  Values  of  n?,  the  nonlinear  index  arising  from  molecu¬ 
lar  reorientation  and  redistribution,  calculated  from  Eq.  (41)  for 
various  liquids.  Values  of  the  threshold  power  for  self-focusing  /’, 
calculated  using  these  values  of  m  in  Eq.  (7),  with  X,  =  6941  A. 
Values  of  rim,  the  nonlinear  index  of  Debye  for  molecular  reorien¬ 
tation  alone,  and  given  by  the  term  in  (41)  lowest  order  in  the 
density.  The  optical  polarizability  data  used  in  calculations  was 
taken  from  Table  142072  o|  Ref.  ft2)  whenever  possible,  and 
otherwise  from  Ref.  (13)  or  as  noted, 


Liquid 

wiX  10** 
esu 

/Vilf)  calc, 
for  6943  A 

ttirXlO1* 

esu 

1  -chlnronaphthalcne  * 

167h 

1.0 

Lead  lrlraisopropvl* 

tit 

1.5 

0 

Tin  tetra  2-methjlbulyl* 

too 

1.5 

0 

Tin  telrapental* 

108 

1.6 

0 

CS, 

96 

1.7 

22 

Tin  tctrabutyl* 

89 

1.9 

0 

SiBr, 

82 

2.0 

0 

Lead  tetraethvl* 

81 

2.1 

0 

Benzoylchloride 

73» 

2.3» 

Nitrobenzene 

69 

2.4 

8.8 

Rromnbenzene 

68' 

2.5 

6.4' 

Acetophenone 

S6>< 

3.1d 

Chlorbenzene 

54 

3.1 

6.3 

Toluene 

45 

3.8 

5.6 

SnCI,i 

44 

3.8 

0 

Benzene 

37 

4.5 

4.4 

Sn(CHj)i 

37 

6.2 

0 

CCI4 

24 

7.1 

0 

Chloroform 

t9 

9.0 

1.3 

Water 

11” 

15b 

Acetone 

6.7 

26 

06 

Acetic  acid 

5.7b 

30*> 

Liquid  CH« 

3.0 

56 

0 

Liquid  Hi 

0.024* 

7000* 

0.010 

•  A  molecule  whose  largest  dimension  ts  estimated  to  be  larger  than  • 
hence,  the  listed  index  and  power  values  are  less  accurate  for  this  than  the 
smaller  molecules. 

k  Measured  values  of  a  not  available ;  m  calculated  using  A  estimated 
from  the  ac  Kerr  effect  data  of  Ref.  (a)  and  the  theory  of  Ref.  (8). 

•  Value  of  A  estimated  by  sealing  values  In  Table  142072  of  Kef.  (12)  to 
those  of  Table  142071. 

4  Estimated  by  assuming  molecule  is  symmetric  (A— 0). 

•  This  value  Is  less  accurate  than  others  because  the  quantum  effects  that 
wc  have  ignored  are  Important  in  liquid  Ek 


index  does  not  contain  any  effects  of  macroscopic 
density  changes  such  as  would  arise  after  a  certain  time 
from  electrostriction  and  heating  in  a  nonuniformly 
irradiated  sample  or  in  a  compressible  sample  uniformly 
irradiated.  The  expression  (41)  contains  only  the  rela¬ 
tively  quickly  established  nonlinearilies  due  to  molec¬ 
ular  reorientation  and  redistribution. 

In  order  to  evaluate  (41),  a  knowledge  of  the  molecu¬ 
lar  volume  Jr  is  required.  This  parameter  also  occurs 
in  the  theories  of  the  fluid  viscosity,  the  heat  conduc¬ 
tivity,  and  virial  coefficients.  An  examination  of  data 
on  representative  fluids  in  each  of  these  areas  gives 
values  for  r  which  vary  typically  by  ±10%  for  a  given 
fluid  and  which  fall  within  ±15%  or  2.5/p.  Since  data 
on  r  are  not  available  for  most  of  the  fluids  of  interest 
to  us,  wc  have  used  the  value  pi  =  2.5  throughout  our 
numerical  evaluations  of  n2.  Fortunately,  n,  is  not  very 
sensitive  to  deviations  of  pr  from  2.5,  and  we  feel  that 
the  errors  arising  from  this  approximation  are  less  than 
from  other  sources. 

The  values  of  u2  that  result  from  using  this  approxi 
mation  and  polarizability  data  appropriate  for  the 
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wavelength  0.589  p  of  the  sodium  D  line  are  given  for 
24  liquids  in  the  first  column  of  Table  I.***1*  The 
threshold  powers  P.  for  self-focusing  estimated  by 
using  our  values  for  n2  in  (7)  are  listed  in  the  second 
column.  Also  given  in  the  last  column  for  comparison 
are  the  values  rtjn  which  would  be  obtained  from  the 
low  density  formula  of  Debye  for  the  nonlinear  index 
[i.e.,  the  first  term  on  the  RJIS  of  (41)]  which  neglects 
molecular  redistribution. 14 

VI.  EXAMINATION  OF  ACCURACY  AND 
DISCUSSION 

The  inequality  (40)  is  equivalent  to  the  statement 
that  the  Clausius-Mossotti  function  («-l)/[p(€+2)] 
is  always  larger  than  4xa/3  for  a  classical  fluid  of 
molecules  of  fixed  linear  polarizability,  regardless  of 
the  intermolecular  forces,  provided  that  the  two- 
particle  correlation  function  depends  only  on  the  inter¬ 
molecular  spacing  and  becomes  constant  for  large 
spacing.  The  fact  that  many  fluids  have  been  found 
which  slightly  disobey  (40)  demonstrates  that  at  least 
one  of  these  assumptions  is  not  entirely  valid.11  How¬ 
ever,  the  observed  deviai.on  of  the  ClausiuS-Mossotti 
function  from  4t«/3  is  rarely  more  than  a  few  percent 
in  liquid.1*  This  suggests  that  our  nonlinear  index 
should  be  accurate  to  within  10  to  50%,  depending  on 
the  molecule.  However,  in  arriving  at  (41),  several 
added  approximations  have  been  made  so  that  this 
absolute  accuracy  may  not  be  attained.  The  result  for 
nt  is  sensitive  to  the  forms  of  the  two-,  three-,  and  four- 
particle  correlation  functions.  The  first  of  these  we 
have  approximated  by  the  step  function  (37),  and  even 
the  best  parameter  r  to  use  in  (37)  is  uncertain. 
Probably  least  certain  of  ali  the  approximations  is  the 
Kirkwood  superposition  approximation  for  the  three- 
and  four-partidc  correlation  functions.  Also,  for  those 
“large”  molecules  listed  in  Table  I  whose  largest  di¬ 
mensions  are  probably  larger  than  extra  errors  are 
introduced  because  such  molecules  are  clearly  not 
representable  as  hard  spheres  with  ideal  polarizable 
dipoles  at  their  centers.  Furthermore,  in  the  listed 
threshold  powers  there  are  the  added  errors  inherent 
in  (7)  which,  when  applied  to  a  given  laser  beam,  could 
easily  err  by  an  order  of  magnitude.  Despite  these  and 
other  difficulties,  we  believe  there  are  some  significant 
correlations  between  our  theoretical  values  of  n j  or  P, 

11  Lendoit-Bomstein,  Zahlenwcrte  und  Function#*  1/3 
(Springer-Verlag,  Berlin,  1962),  p.  509. 

»  Dictionary  of  Orranic  Compounds  (Oxford  Universily  Press, 
New  York,  1965). 

>‘Y.  R.  Shen,  Phys.  Letters  20,  378  (1966),  also  uses  this 
formula  for  nm  lo  tabulate  nonlinear  indices.  However,  his  values 
often  differ  from  ours  by  a  factor  ~2  because  he  employs  the  dc 
Kerr  coefficients,  with  data  on  permanent  dipole  moments  and 
the  theory  of  C.  V.  Raman  and  K.  S.  Krishnan,  Phil.  Mag.  3,  724 
(1927),  to  determine  6,  rather  than  the  optical  values  measured 
directly  by  Rayleigh  scattering.  This  latter  theory  was  shown  by 
Raman  and  Krishnan  to  be  commonly  in  error  by  2  or  more. 

•*  C.  M.  Knobler,  C.  P.  Abiss,  and  C.  J.  Pings,  J.  Chem.  Phys. 
40,  2200  (1964). 


and  the  available  experimental  results  on  self-focusing. 
The  reasons  for  this  and  some  detailed  comparisons  with 
experiment  are  given  below. 

From  recent  time-resolved  studies  of  the  formation 
of  self-t rapped  filaments  of  an  unfocused  ruby  laser 
beam  in  CSs  and  nitrobenzene,  it  was  found  that  the 
filaments  do  not  persist  for  more  than  a  few  nano¬ 
seconds  even  when  the  laser  pulse  is  several  tens  of 
nanoseconds  long.1*  7n  all  of  the  experimental  studies 
of  self-trapping  to  date,  unfocused  beams  of  cross 
section  of  the  order  of  1  mm  or  larger  have  been  used. 
Therefore,  during  the  lifetime  of  a  “filament  ’’  a  sound 
wave  could  travel  only  a  small  fraction  (<10~*)  of  the 
width  of  the  parent  beam.  This  suggests  that  macro¬ 
scopic  density  changes  and  the  effects  (such  as  elect  re¬ 
striction)  which  produce  them  may  not  have  time  to 
affect  the  index  in  many,  if  not  all,  materials  studied 
to  date.  Further  evidence  against  density  change  effects 
comes  from  the  observed  instability  in  the  nonlinear 
propagation  of  a  circularly  polarized  beam,  an  insta¬ 
bility  which  would  not  exist  if  density  changes  alone 
altered  the  index.'*  It  is  true  that  stimulated  Brillouin 
scattering  may  accompany  or  precede  self-focusing,  but 
the  usual  theory  suggests  that  the  density  changes 
associated  with  this  are  not  such  as  to  case  an  index 
change  at  the  parent  frequency,  although  they  may  of 
course  contribute  a  loss  which  can  effect  self-focuring. 
Therefore,  it  seems  reasonable  to  expect  that  a  com¬ 
parison  between  present  data  and  our  theory  which 
omits  effec.  of  macroscopic  density  changes  will  j  .'eld 
some  information  on  the  accuracy  of  the  theory  and  of 
the  conclusions  that  it  suggests.  Such  a  comparison 
and  suggestions  for  more  definitive  experimental  checks 
of  molecular  redistribution  effects  are  now  given. 

In  the  first  data  on  self-trapping  of  Hauchecorne  and 
Mayer,*  nitrobenzene  was  found  to  self-focus  about  the 
same  fraction  of  an  i  cident  beam  as  does  1-chloro- 
naphthalene.  This  suggests  that  the  threshold  powers 
are  about  the  same  in  these  two  materials,  as  has  been 
verified  more  recently  by  Bret  el  al  •  Therefore,  the  ratio 
of  the  threshold  for  1-chloronaphthalene  to  that  of 
nitrobenzene  inferred  from  Table  I  is  two  times  tor 
low.  However,  1-chloronaphthalene  is  a  “large”  mole¬ 
cule  (its  maximum  dimension  is  greater  than  p~in)  for 
which  our  model  of  a  point  i  role  at  the  center  of  a 
spherical  molecule  is  obviously  bod.  This  suggests  that 
the  critical  powers  listed  for  the  other  “large”  molecules 
indicated  in  Table  I  may  also  be  too  low  relative  to  the 
other  powers  listed. 

Shen  and  Shaham*  fo^nd  that  the  threshold  powers 
for  self-focusing  in  nitrobenzene,  acetophenone,  ben¬ 
zene,  and  water  increased  in  that  order.  This  is  in 
agreement  with  the  order  of  Table  I. 

11 D.  H.  Close,  C.  R.  Guiliano,  R.  W.  Heliwxrth,  L.  D  Hess, 
F.  J.  McClung,  and  W.  G.  Wagner,  IEEE  J.  Quant.  Electron. 
QE-2,  No.  9,  553  (1966). 
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’  illt  riMid  am!  Blovmbvtgvn  found  that  the  threshold 

.1  its  of  US,,  nitrobenzene,  bromoben/.enc,  and 
aeetone  increased  in  that  order,  in  agreement  with 
Table  1-  They  also  failed  to  observe  self-trapping  in 
water  and  l’l. "I*.  However,  from  a  more  detailed 
eviniit’.ttinn  of  Raman  scattc  ing  near  threshold, 
Bret  ( l  al.  have  found  evidence  for  self-trapping  in  CC'I4 
at  powers  w  here  there  was  none  for  acetone,  chloroform, 
and  acetic  acid.’  Table  I  is  consistent  with  the  latter 
possibility. 

Other  measurements  of  ret  et  al.  give  the  threshold 
powers  Coproduee  stimulated  Raman  scattering  in  CSj, 
nitrobenzene,  benzoylchloride,  1-chloronaphthalcne, 
toluene,  benzene,  and  C'C'l4  to  be  in  the  ratios  1: 2:3:3: 
6:10:300*  Since  Bret  et  al.  present  indirect  evidence 
that  self  trapping  occurs  in  each  of  these  liquids,  we 
tnav  suppose  that  this  ordering  is  also  that  for  the  self- 
t  rapping  threshold  powers.  Again,  except  for  the  “large" 
molecule  1  cbloronaphthalene,  the  ordering  is  consistent 
with  that  in  Table  I.  The  connection  between  these 
threshold  powers  covering  a  power  range  of  300  and 
the  self-trapping  thresholds  for  beams  of  a  fixed  spatial 
pattern  is  not  clear  enough  to  warrant  studying  a 
numerical  comparison,  especially  since  the  experimental 
procedure  was  not  given. 

From  a  study  of  ‘.he  cell  length  dependence  of  stimu¬ 
lated  Raman  thresholds,  Wang  has  inferred  that  the 
trapping  threshold  powers  for  nitrohenezene,  toluene 
and  benezene  arc  in  the  ratios  1:2.9:3.4.  This  may  be 
compared  with  corresponding  ratios  1:1.6: 1.9  from 
Table  I.  We  do  not  feel  that  this  difference  represents 
a  fundamental  disagreement,  but  is  probably  repre¬ 
sentative  of  the  experimental  and  theoretical  errors  in 
such  ratios. 

If  even  rough  comparisons  are  to  be  made  between 
the  absolute  values  of  the  threshold  powers  for  self- 
focusing  predicted  by  (7)  and  obser/ation,  a  highly 
degenerate  diffraction  limited  beam  of  smooth  wave- 
front  must  be  employed.  Otherwise  corrections  tc  (7) 
discussed  by  Wang4  must  be  devised.  The  only  pub¬ 
lished  efforts  to  st  udy  self-focusing  w ith  a  beam  approxi¬ 
mating  the  ideal  conditions  are  those  of  Garni  ire  et  al.e 
Thev  found  beam  constriction  to  set  in  at  25 ±5  kW 
in  nitrobenzene,  an  order  of  magnitude  higher  power 
than  given  b'  our  value  of  n2  used  in  (7).  With  a  multi¬ 
mode,  multi-lobe  beam,  Wang  inferred  from  Raman 
data  that  the  self-reusing  in  nitrobenzene  set  in  at 
19  kW,4  again  a  much  higher  power  than  in  Table  I. 
Such  discrepancies  could  not  be  said  to  be  surprising 
in  view  of  the  combined  absolute  errors  of  (7)  and  (11). 
However,  there  is  some  evidence  that  it  is  premature 
to  a  -sume  that  this  discrepancy  is  entirely  established, 
first,  Chino  and  Gnrmirc  have  reported  that  very  small 
filaments  have  been  discovered  subsequently  to  be 
present  while  the  overall  constriction  of  the  beam,  on 
which  the  previous  threshold  estimate  was  based,  is 


just  beginning.17  'mis  sub  structure  in  the  beam  sets 
in  at  an  us  yet  undetermined  lower  power  level  than  the 
large  scale  beam  constriction.  Secondly,  Emmett  has 
observed  in  nitrobenzene  what  appear  to  be  self -trapped 
filaments  in  the  stimulated  Raman  scattering  at  90° 
from  a  beam  that  has  been  focused  onto  a  perpendicular 
line."  These  “filaments”  vary  in  number  with  the  beam 
power  as  do  the  usual  parallel  filaments.  However, 
single  “filament"  outputs  of  between  5  and  10  kW 
have  been  observed,  in  rough  agreement  with  the  values 
of  Table  I."  Evidently,  further  work  is  required  to 
establish  definitive  estimates  of  thresholds  for  the 
self-trapping  of  highly  degenerate  diffraction  limited 
beams. 

There  are  two  kinds  of  measurements  which  have 
been  performed  that  do  not  involve  «2  or  self-focusing, 
but  which  may  give  information  about  redistribution 
effects,  first,  Mayer  and  '’ires  have  measured  the 
optical  birefringence  induced  in  one  beam  traversing  a 
material  fc  •  another  strong  linearly  polarized  congruent 
beam  of  a  different  wavelength.*'1*  They  describe  their 
results  by  an  optical  Kerr  constant  B0,  which,  if  the 
low  density  theory*  were  correct,  would  equal  3h2/>/ 
(2 X i ) ,  vhvre  At  is  the  free  space  wavelength  of  the 
birefringent  beam.50  However,  when  molecular  redis- 
tributiot.  effects  at  liquid  densities  are  also  considered, 
a  diffemu  relation  between  B0  and  «2  will  result.  Un¬ 
fortunate!,  our  theory  of  «2  is  not  capable  of  giving 

“  “twv  beam”  property  Bo,  but  it  could  be  extended 
to  do  so.  After  such  an  extension,  a  comparison  of 
predicted  values  of  B0  with  those  that  have  been  meas¬ 
ured  would  probably  make  a  much  more  accurate  test 
of  theory  than  comparisons  of  n2,  because  the  measure¬ 
ments  of  Bt  are  probably  more  accurate  than  those 
of  n2.*1 

Secondly,  Maker  et  al.  have  measured  the  change  in 
the  state  of  polarization  of  a  strong,  elliptically  polar¬ 
ized,  b< am  as  a  function  of  beam  power  and  distance 


11 R.  Y.  U'hiao  and  F..  Garmire  IEEE  J.  Quant  Electron. 
QE-2,  No.  9,  467  (t966). 

"  J.  I..  Emmett  (private  communicalion). 

’*  G.  Mayer  and  F.  Gires,  Compt.  Rend.  258,  2059  (1964). 

*  A  comparison  of  the  observa'ions  of  Mayer  a,,d  Gires  with 
predictions  of  the  low -density  theory  of  Ref.  8  shows  general  dis¬ 
agreement.  Reference  9  reports  that  /J»Xt0*  esu  is  observ'd  to  be 
42,  29,  9.5,  4.0,  0.5,  t.6,  and  0.4  for  the  liquids  CSi,  nitrobenzene, 
toluene,  benzene,  CCb,  cholofoim,  and  cyclohexane,  respectively 
(at  Xi  =  488  mri).  Using  the  optical  polarizabilities  obtained  from 
Rayleigh  scattering  measurements  (Ref.  t2)  in  the  low  density 
theory  (Ref.  8),  one  would  predict  69,  27,  18,  14,  0,  4.0,  and  t.5 
for  these  values. 

51  This  extension  of  theory  would  also  give  predictions  for  the 
static  (dc)  Kerr  constant  B.  for  nonpolar  molecules  which  could 
be  compared  with  a  wealth  of  xisting  data.  These  data  are  in 
w  ide  disagreement  with  the  predictions  of  the  classical  low  density 
theory  4rU:.  8).  For  example,  the  measured  values  of  if.XtO’ esu 
for  US;,  benzene,  CCIi,  liquid  Nr,  cyclohexane,  and  liquid  Hj,  at 
Xi  =  546  m*i  a:e  55. 5,  4 . t ,  0.84,  0.8,  0.59,  and  0.54,  respectively 
[front  fandolt-Uomstein,  Zahlenu-trlt  nnd  lunktionen  U/S 
(Sprincer-Yerlag,  Berlin,  1962),  Chap.  5,  pp.  849-855],  Using 
the  {Hiiarizability  data  of  Kef.  12  in  the  low  densitv  theory  (Kcl. 
8),  one  would  predict  Ml,  t2,  «,  11.4,  t  4,  and  11.027,  rcjeclivelv, 
for  tluse  values. 


K 


111  1.1.  VV  A  K  III 


164 


ilniu;  i ho  Im.iiii.™  Hum1  incaviirtiiicnls  luvc  been 
lciYiuly  refuted  1  tv  Wang  and  Kacette.-1  Their  result'! 
um‘  given  in  terms  of  a  con>tant  li  (real  in  the  absence 
of  loss5  which  would  ei|iial  3;j;i_.  cir  if  the  low  density 
theory '  which  mgleits  molecular  reilistriliution  were 
ipplicable.  Inti  w ill  be  relatn!  Ii  i,  ■  d’lTerently  at  liquid 
densities  where  ninleculai  udis! dilution  is  important. 
\gain,  when  our  theory  is  extended  to  predict  It  values, 
comparison  w  it  li  these  experiments  w  ill  be  illuminating. 

Titer'.-  is  a  measurement  of  ;i5  which  could  be  per¬ 
formed,  but  which  does  not  involve  self  focusing  and 
hence  might  he  more  accurate.  The  change  in  index  for 
a  weak,  plane  polarized,  beam  congruent  with  a  strong 
beam  of  the  same  polarization  hut  of  a  different  wave¬ 
length  would  equal  n:E:,  where  E-  is  tile  mean  square 
amplitude  of  the  strong  beam,  provided  that  the  two 
wavelengths  are  not  too  widely  separated.  This  change 
in  index  could  be  measured,  for  example,  by  putting 
the  strongly  irradiated  material  in  one  amt  of  a  Mach- 
Zehnder  interferometer,  and  using  the  weak  beam  at  a 
different  wavelength  as  the  interferometric  source. 

Stimulated  llrillouin  scattering,  turbulence,  impuri¬ 
ties,  bubbles,  and  molecular  relaxation  processes,  all 
of  which  tend  to  prevent  the  liquid  from  reaching  an 
equilibrium  state  in  the  strong  fields,  may  affect  experi¬ 
mental  results.  Elucidation  of  these  phene  ;na  will 
doubtless  follow  a  better  understanding  of  the  equi¬ 
librium  state  itself,  but  we  will  not  attempt  to  analyze 
them  here. 

In  conclusion,  since  there  are  no  general  relations  for 
liquids  among  the  coefficients  n j,  Bo,  and  B,  information 
concerning  the  validity  of  approximations  such  as  the 
Kirkwood  approximation  in  dealing  with  molecular 
redistribution  effects  can  be  obtained  as  yet  only  from 
comparing  the  present  theory  with  experiments  which 
measure  Hi-  When  the  theory  is  extended  to  treat  Bo 
and  B,  additional  checks  of  the  approximations  will 
become  ax’ailable.  However,  with  only  the  relatively 
crude  comparisons  between  the  available  results  on 
self-trapping  and  the  theory  that  we  lmve  made  aboxe, 
we  arc  led  to  conclude  that  molecular  redistribution 
effects  are  generally  important,  and  no  d<  ubt  dominate 
the  nonlinear  index  of  symmetric-molecule  liquids. 
Observations  of  lf-t rapping  in  such  symmetric- 
molecule  liquids  as  SiBr«  and  :nCl«  will  provide  im¬ 
portant  checks  of  the  theory  of  the  nonlinear  index. 
Further  efforts  to  obtain  smooth,  diffraction  limited, 
highly'  degenerate  beams  would  aid  in  establishing  the 
absolute  magnitude  of  the  nonlinear  index  ut  in  various 
materials.  The  measurement  of  the  index  change  for  a 
weak,  linearly  polarized,  beam  congruent  within  a 
material  to  a  parallel -polarized  strong  beam  of  different 
wavelength  may  yield  an  even  more  accurate  value  of 

for  that  material. 


”  P.  1).  Maker,  R.  \X'.  Trrhune,  and  C.  M.  Savage,  Pliys.  Rev. 
I  .oilers  12.  507  (1964) 

11 C.  C  Wang  and  (.«  W.  Race-lie,  J  (Juan 1 .  1'JrClron.  OP-2, 
No  4  54  (t%6). 


APPENDIX 

llcie  we  outline  the  calculation  of  U77  =  <(«2J)) 
-  «rj)/.  We  will  use  the  form  (.57)  for  the  probability 
\'~-p(ti  —  TiVVit/Vj  that  if  a  certain  molecule  is 
in  rfVi,  another  specified  molecule  is  in  <Pr7.  We 
will  use  Kirkwood’s  superposition  approximation  for 
I  -3/»a(r,,r2.ra)  =  />aC123)  [and  for  V~*pt( ri,r,,fi,r4) 
=  /b(1234)]  which  is  the  probability  (per  unit  volume 
for  each  molecule)  for  finding  specified  molecules  at 
tn  tji  (and  r<)  simultaneously : 

/*,(123)«/»(12)/»(23)*(13) ,  (Al) 

^(1234)«/>(12)/>(i3)/K14)/K23)M24)£(34).  (A2) 

Here  p(\2)  is  an  abbreviation  for  p( r»),  and  we  have 
ignored  terms  of  order  A’_:  smaller  than  the  leading 
terms  of  (Al)  and  (A2).  (Terms  of  order  Ar  must  be 
kept,  as  the  terms  in  l\7  which  arc  proportional  to  A‘* 
w  ill  cancel.) 

The  integrals  in  4«h,j))  fall  naturally  into  four 
gioups,  which  we  call  Ii,  I  j,  /  j,  and  /«,  depending  on 
whether  1-,  2-,  3-,  or  4-particle  correlations  are  required : 

4  ««:’)>=£/.■  (A3) 

1-1 

where 

/,=£„  «(E'-a<‘-D--«»  E,)5j)=0  (A4) 

/,=  2i;,.«(E'-«-D-'-a'-E7»  (A5) 

fi  =  4 ((£'•«’  ’  D"' «'■  E’E'' «'•  !)’*■  a’- E'))  (A6) 

and 

<A7) 

Here  the  summed  indices  are  alway  s  unequal  in  any  / 1. 
The  1 1  is  zero,  of  course,  because  D‘“‘  is  defined  to  be 
zero.  The  angular  averages  in  each  of  the  remaining 
integrals  are  easily  performed  and  give 

I}  =  2(.tI■.',  Enr  [(I+A,'4)l<(/2„-)5) 

+  (A/2+  llAV16)CI‘rD'"-D")]  (A8) 

4o4/i'‘(l  +  A)  (I)..“’I),/-<)  (A9) 

and 

/4=«<F'«  {D.S’D"*').  (AIO) 

Here  “  IT”  denotes  the  "trace”  or  sum  over  diagonal 
space  elements,  and  llie  single  bracket  (  )  denotes 

the  average  over  spatial  configurations. 

Using  the  form  (37)  for  the  two  particle  distribution 
function  and  (20)  for  D"  gives  (neglecting  ternis~A'°) 

Z..((P..I‘')7)=X1V  ‘Hr/’  sin  W*\  r'dr 
J  0  j  d 


X  (1  — 3  co^VOV  E  =  .Vp64jr/(45r).  (All) 
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Similarly, 

v;„.  <  IVD-  D-')  XiV^'l2rf  sin ^#  /  r1  dr 

J  o  J  J 

X  ( 1  +  3  cosV)r~* = A’p-liir3/  (4r) .  (A12) 

Therefore, 

/ .  =  A  V/-;  V4)r  (52  + 1 46A+ 167  A3)/  (45pr) .  (A  14) 

To  perform  the  average  in  (A9)  to  order  A',  we  write 

/>(r)  =  (1  —/;(/))  (1  —  r/ 1') ,  (A  14) 

where  A(r)=l  for  r<d  and  is  zero  otherwise.  Then, 
with  (Al),  we  obtain  to  order  A’ 


The  right  hand  side  separates  naturally  into  terms 
which  we  call  K,  where  »  =  (>,  1,  •  •  •,  4  is  the  number  of 
explicit  It  factors  in  the  integrand  of  the  term.  The  A"o 
term  is  integrated  trivially  and  has  a  part  of  order  A’2 
which  cancels  the  A'3  part  of  ((«*))*,  and  a  part  of  order 
A-,  which  combines  wi  h  ((ih))1  from  (41)  (corrected 
by  the  now  important  factor  (1  —  r/l')_l)  to  give 

A'o-4((«2))2=  -4AVf3(l-pr).  (A  19) 

There  are  four  terms  in  A'i  (there  being  four  possible 
terms  linear  in  one  of  the  h’ s)  each  of  which  integrates 
trivially  to  —Ftp1  giving 


-  A’3 F  j  (Pr/Pr-DitDail— /»i; ) ,  (A15) 

where  Du^l)l,',{l—ltit),  denotes  //(r2-r3),  etc. 
The  integrals  here  are  easily  performed  using  the 
theorem  that  the  integral  J  of  /),»'*  over  any  volume  Q 
is  related  to  an  integr  !  over  the  surface  S  bounding 
Q  by 

J=  [  [  J*riti  itl..-irtrt.  (A  16) 

J  Q  J  S 

(If  S  is  ellipsoidal  and  r,  is  inside  it,  recall  that  J=L. 
If  5  is  spherical  and  ri  is  outside  it,  then  J-QD 
where  u  is  the  center  of  the  sphere  and  Q  is  the  nu¬ 
merical  value  of  the  volume  bounded  by  S.)  The  part 
of  (A15)  independent  of  itn  is  obviously  A’p2/2.  With 
the  aid  of  (A16),  the  part  of  (A  15)  involving  It u 
integrates  to  A'ff2ry3  whence 

/i-4AV£'V(l+A)(P+2«3/3) .  (A  17) 

To  evaluate  lt  we  use  (A2)  to  obtain 

£„.i,  (D.S’D,,3') 

=  ,V(A’—  1)  (A*-  2)  (A*  -  4)  1-J(1-  r/  V)~* 

X  (1  —  /r  i  j)  (1  —  hi*) 

=  £A\.  (A  18) 

1-0 


Ki=  —4Fp*tN.  (A  20) 

There  are  six  ways  to  form  integrands  quadratic  in  the 
It's  in  (A18),  four  of  which  have  one  space  point  as  an 
argument  of  both  It’s  and  are  equivalent  to  each  other. 
The  remaining  two  have  no  space  points  as  arguments 
of  both  It  factors  and  are  equivalent  to  each  other. 
Each  of  the  first  four  terms  is  seen  to  be  iderti. — !v 
ze;o  with  the  aid  of  (A16),  and  the  last  two  can  be 
straightforwardly  integrated  with  (A16)  to  give 

A‘j  =  0.5066irp!rA’ .  (A21) 

Almost  the  same  straightforward  integrations  are 
required  for  each  of  the  four  equivalent  terms  which  are 
cubic  in  the  It’s,  and  one  finds  that 

A’j=-0.11713jrp3nV.  (A22) 

We  have  not  been  able  to  integrate  the  remaining 
term  Kt  analytically.  We  estimate  that  it  is  an  order 
of  magnitude  smaller  than  A'j  because  the  four  It  factors 
combine  to  l:~>it  the  volume  accessible  to  the  points  of 
integration  much  more  than  in  A"Y  Since  ATj  is  already 
a  minor  contribution  to  the  final  result,  we  feel  well 
within  the  other  limits  of  our  calculation  to  neglect 
A'«  in  the  final  result: 

It— 4((/ij))5~  Ap3(  —  4F-f-3.84pr) .  (A  23) 

Combining  all  the  above  terms  gives  the  l/jj  of  (38). 


L 

l 

L 


* 


. . . i-Tn—n - - iiiiiiiiiiiiiiiiiiiiiiniiiiniiiimiimiiiiiinitiiiiiiiilitiiiirillliiininiitrilliilliiitinntitiiiiiirntiiiiiiiiiiiiiiiiiiiiiiitiiiiiii  - - - m - r~ - rr - - - — - - - 


SSI  FI  ED 
Security  Classification 


DOCUMENT  CONTROL  DATA  -  RAD 


?b  CROUP 


(Saairtty  ctaaallicatton  ol  till*  body  ol  abatraet  and  inaaiin|  annotation  null  ba  antarad  whan  tha  oratat'  tapott  ia  ataaallladj 


I  ORIGINATING  ACTIVITY  (Corporata  author)  la  REPORT  lECuRi  ty  c  laisipication 

Hughes  Research  Laboratories  Unclassified 

3011  Malibu  Canyon  Road 
Malibu,  California  90265 


i  report  title 

Study  of  Fundamental  Requirements  for  High  Brightness  Raman 
Laser  Sources 


4  DESCRIPTIVE  NOTES  (Typ*  ot  rtporf  #nd  inc/utiv*  d#f##J 

Semiannual  Technical  Summary  Report,  1  October  1966  -  31  March  1967 


5  AUTNORfSJ  (Lm»t  nm »•.  tint  nmmm.  Initial ) 

McClung,  F.  J. 


«.  REPORT  DATE 


May  1  967 


•  a.  CONTRACT  OR  GRANT  NO. 

NO  0014-67 -C-0206 

b  PROJECT  NO. 


7a  TOTAL  NO  OP  PAGES  I  7  6.  NO.  OP  REPS 


■  a.  ORIGINATOR'S  REPORT  NUMBERfSj 


ARPA  Order  No.  306 


16.  OTHER  REPORT  NOffJ  ( A  ny  othar  nunbcni  A«l  (My  6* 
Ell#  report? 


10  A  VA  IL  ABILITY/LIMITATION  NOTICES 


Qualified  requesters  may  obtain  copies  of  this  report  from  DDC 


II.  SUPPLEMENTARY  NOTES 


13  ABSTRACT  j- 


12.  SPONSORING  MILITARY  ACTIVITY 

Physical  Sciences  Division 
Office  of  Naval  Research 


IVFM  illliu  w»f 


We  have  studied  the  SRS  in  liquids.  The  production  of  high  brightness  Raman 
devices  is  possible  in  liquids;  however,  the  brightness  gains  over  existing 
laser  devices  does  not  appear  to  be  very  large.  We  must  exercise  care  in 
selection  of  the  material  to  achieve  any  semblance  of  gain. 

The  case  for  gases  is  much  better.  Many  unexplained  processes  occur  in 
H2  and  possibly  other  gases  as  well.  The  brightness  gains  that  could  be 
produced  have  not  yet  been  determined,  and  the  main  effort  is  now  con¬ 
centrated  on  this  problem.  We  are  also  investigating  the  strange  behavior 
of  in  the  usual  SRS  experiments. 


DD 


FORM 

1  JAN  64 


1473 


UNCLASSIFIED _ 

Security  Classification 


UNCLASSIFIED 


Security  Classification 


La  s  c  r  s 

Raman  Scattering 
Nonlinear  Optics 
Liquids 
Gase  s 

H? 

Self  Focusing 


INSTRUCTIONS 


I.  ORIGINATING  ACTIVITY:  Enter  the  name  and  addreaa 
of  the  contractor,  aubcontractor,  grantee.  Department  of  De- 
fenae  activity  or  other  organisation  ( corporate  author)  iaaulng 
the  report, 

2a.  REPORT  SECURITY  CLASSIFICATION:  Enter  the  over¬ 
all  aecurity  classified  on  of  the  report.  Indicate  whether 
"K.atricted  Data"  ir  include*!  Marking  la  to  be  in  accord¬ 
ance  with  appropriate  aecurity  regulations. 

2b.  GROUP:  Automatic  downgrading  ia  apecified  in  DoD  Di¬ 
rective  5200. 10  and  Armed  Forcea  Industrial  Manual.  Enter 
the  group  number.  A!ao,  when  applicable,  ahow  that  optional 
markinga  have  been  used  for  Group  3  and  Group  4  as  author¬ 
ised. 

3.  REPORT  TITLE:  Enter  the  complete  report  title  in  all 
capital  letter  a.  Titlea  in  all  caaea  should  be  unclassified. 

If  a  meaningful  title  cannot  be  aelected  without  classifica¬ 
tion,  ahow  title  classification  in  all  capitals  in  parentheaia 
immediately  f<  .wing  the  title. 

4.  DESCRIPTIVE  NOTES:  If  appropriate,  enter  the  type  of 
report,  e.g,,  interim,  progress,  summary,  annual,  or  final. 

Give  the  inclusive  dates  when  i  specific  reporting  period  is 
covered. 

5.  AUTHOR(S):  Enter  the  name(a)  of  author!*)  »a  shown  on 
or  in  the  report.  Enter  last  name,  first  name,  middle  initial. 

If  military,  show  rank  and  branch  of  service.  The  name  of 
the  principal  author  is  an  absolute  minimum  requirement 

6.  REPORT  DATE:  Enter  the  date  of  the  report  aa  day, 
month,  year,  or  month,  year.  If  more  than  one  date  appears 
on  the  report,  uae  date  of  publication. 

7a.  TOTAL  NUMBER  OF  PAGES.  The  total  page  count 
should  follow  normal  pagination  procedures,  Le.,  enter  the 
number  of  pages  containing  information. 

7b.  NUMBER  OF  REFERENCES:  Enter  the  total  number  of 
references  cited  in  the  report. 

8a.  CONTRACT  OR  GRAN  f  NUMBER:  If  appropriate,  enter 
the  applicable  number  of  the  contract  or  grant  under  which 
the  report  was  written. 

8b,  8c,  b  8«f.  PROJECT  NUMBER:  Enter  the  appropriate 
military  department  identification,  such  aa  project  number, 
subproject  number,  system  numbers,  task  number,  etc. 

9a.  ORIGINATOR'S  REPORT  NUMBER(S):  Enter  the  offi¬ 
cial  report  number  by  which  the  document  will  be  identified 
and  controlled  by  the  originating  activity.  This  number  must 
be  uniqua  to  this  report. 

9b.  OTHER  REPORT  NUMBER(S).  If  She  report  has  been 
assigned  any  other  report  numbers  (either  by  the  originator 
or  by  the  sponsor),  also  enter  this  number! s). 

10.  AVAIL  ABILITY/ LIMITATION  NOTICES:  Enter  any  lim¬ 
itations  on  further  dissemination  of  the  report,  other  than  those 


imposed  by  security  classification,  using  standard  statements 
such  aa: 

(1)  "Qualified  requesters  may  obtain  copies  of  this 
report  from  DDC." 

(2)  "Foreign  announcement  and  dissemination  of  this 
report  by  DDC  is  not  authorized." 

(3)  "U.  S.  Government  agencies  may  obtain  copies  of 
this  report  directly  from  DDC.  Other  qualified  DDC 
users  shall  request  through 


(4)  "U.  S.  military  agencies  may  obtain  copies  of  this 

report  directly  from  DDC  Other  qualified  users 
shall  request  through 


(5)  “All  distribution  of  this  report  ia  controlletL  Qual¬ 
ified  DDC  users  shall  request  through 


If  the  report  haa  been  furnished  tc  the  Office  of  Technical 
Services,  Department  of  Commerce,  for  sale  to  the  public,  indi¬ 
cate  this  fact  and  enter  the  price,  if  knowiv 

1L  SUPPLEMENTARY  NOTES:  Use  for  additional  explana¬ 
tory  notes. 

12.  SPONSORING  MILITARY  ACTIVITY:  Enter  the  name  of 
the  depert menial  project  office  or  laboratory  sponsoring  (pay¬ 
ing  for)  the  research  and  development.  Include  address. 

13.  ABSTRACT:  Enter  an  abstract  giving  a  brief  and  factual 
summary  of  the  document  indicative  of  the  report,  even  though 
it  may  also  appear  elsewhere  in  the  body  of  the  technical  re¬ 
port.  If  additionsl  space  is  required,  a  continuation  sheet  shall 
be  attached. 

It  is  highly  deairsble  that  the  abstract  of  classified  reports 
be  unclassified.  Each  paragraph  of  the  abstract  shall  end  with 
an  indication  of  the  military  security  classification  of  the  in¬ 
formation  in  the  paragraph,  represented  as  (TS).  (S).  (C).  or  (V). 

There  is  no  limitation  cn  the  length  of  the  abstract.  How¬ 
ever,  the  suggested  length  is  from  ISO  to  225  woids. 

14.  KEY  WORDS:  Key  words  are  technically  meaningful  terms 
or  short  phrases  that  characterize  a  report  and  may  be  used  as 
index  entries  for  cataloging  the  report.  Key  words  must  be 
selected  so  that  no  security  classification  is  required.  Identi¬ 
fiers.  such  as  equipment  model  designation,  trade  name,  military 
project  code  name,  geographic  location,  may  be  uaed  as  key 
words  but  will  be  followed  by  an  indication  of  technical  con¬ 
text.  The  assignment  of  links,  rules,  and  weights  is  optional 


OTO  SM-S51 


UNCLASSIFIED 
Securtty  Classification 


